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Abstract. We use contemporary evolutionary models for Very Massive Stars (VMS) 
to assess whether the Eddington limit constrains the upper stellar mass limit. We also 
consider the interplay between mass and age for the wind properties and spectral mor- 
phology of VMS, with reference to the recently modified classification scheme for 02- 
3.5If*/WN stars. Finally, the death of VMS in the local universe is considered in the 
context of pair instability supernovae. 



1. Eddington limit 

Empirical determination s for the upper stellar mass limit, M max , have led to the adoption 



of M max ~ 150M o (e.g. lFigeril2005l ; IOev & Clarkdl2005h . This limit closely coincides 



with the 'first approximation' for M max ~ 200M o , i.e. the intersection between the 
Eddi ngton limit an d the L oc M 3 mass-luminosity relation for main sequence stars 



(e.g. iMaederl [20091 his Fig 3.6). Therefore, two obstacles need to be overcome for 



these find ings to be r ec onciled with th e claimed 320 M Q initial mass for R136al by 
Crowther et al.l ( 20101) . iMasseyl ( 2011 ) has provided convincing arguments regarding 



the empirical estimates of M max , so here we focus attention on M max . 

Main sequence models for very massive stars (VMS) h ave been calculated u sing 
the Geneva evolutionary code by R. Hirschi and N. Yusof (seelCrowther et all2010l) and 
the Bonn evolutionary code by K. Friedrich (see Gra fener et al. EoT l]). In Figured] we 



present the mass-luminosity rel ation for non-rotating, so lar metallicity ZAMS, incor- 
porating 9-85 M Q models from lMeynet & Maederl (|2000|) . Although x = 3 represents a 
sensible average for L oc M x across all stellar masses, x ~ 2.5 for 10-20 M Q and flattens 
further at higher masses, reaching x ~ 1.5 close to 200 M Q . 
The Eddington parameter, T e , can be expressed as 

L/L Q 



T e = gelg = 3 X 10" 



M/M e 



where q - 0.86 for main sequence hot stars. A decrease in the slope of the mass- 
luminosity relation at very high masses reduces Y e , and so raises ^ max - Of course, Y e 
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Figure 1 . Mass-luminosity relation for no n-rotating, solar metalli city ZAMS stars 
from lMevnet & Maederl d2000l < 85M ) and Crowthe rlstail d2010l > 12OM ). 



increases once a star evolves away from the ZAMS (L/M increases), and since M oc 
L(r g /(1 - T P )) 2 ^ for radiatively driven winds with a CAK po wer index of a ~ 0.6 
dQwockill2003h . str onger winds are anticipated both qualitatively dSmith & Contill2008 ) 
and quantitatively (iGrafener & Hamannll2008h with age. 

Figure [2] compares Y e for ZAMS spanning 10-500 M Q at solar composition. This 
illustrates that the Eddington limit is not approached, so » 500M o . In fact, 

x — > 1 as M — > oo, so the Eddington limit might never be reached for ZAMS stars. 
Once on the main sequence, T e , and in turn mass-loss rates, increase with both age 
and mass, so strong wind signatures may correspond either to a relatively evolved high 
mass star or an unevolved very high mass star. Tabled] compares the influence of mass 
and a ge upon spectral type for the case of the coeval cluster R136a whose age is -1.5 
Myr fcrowther et alj|2010h . 



Table 1 . Influence of mass (vertical) and age (horizontal) upon spectral type for 
the young LMC star cluster R136a 



Initial 


Sp Type 


Sp Type 


Example 


Mass (M Q ) 


(ZAMS) 


(1.5 Myr) 




240 


02 If*? 


WN5h 


R136a2 


140 


02 III? 


02 If* 


R136a5 


100 


02-3 V? 


03 III(f*) 


R136a7 


50 


03 Vz? 


03V 


[HSH95] 50* 



$: iHunter et all dl995L HSH95) 
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Figur e 2. Eddington p arameter, T e , for no n-rotating, solar metallicity ZAMS stars 
from lMeynet & Maederi (l2000l < 85M Q ) and lCrowther et all (I20101 > 12OM ). 



2. Transition Of/WN stars 



From the previous section, high r e 's develop either in young, very high mass stars or 
evolved lower mass stars. Spectroscopic signatures of strong winds in early type stars 
include Hen A46S6 and/or Ha emission, corresponding to OBA supergiants or Wolf- 
Rayet stars in the case of very strong emission features. From Table [TJ the current 
spectral type of R136a2 is WN5h, while its ZAMS spectral type may have resembled 
an 02 supergiant. Had we witnessed R136 perhaps 0.5 million years ago, it would 
have exhibited an intermediate spectral type. Indeed, a hybrid 03 If*/WN category - 
sp ectroscopically intermediate between early O stars and WN stars - was introduced 
bv lWalbornl(ll982h . 

Following the extension of the MK sequence to 02 dWalborn et al.ll2002h and re- 
vision s to WN classifications, this has been refined recently by ICrowther & Walbornl 
(1201 lh . Spectroscopically the morphology of H/3 is ke y to 02-3.5 If*, 02-3 . 5 If*/W N 
or WN subtypes, while a qualitative interpretation led ICrowther & Walbornl ( 2011 ) to 
conclude that most 02-3.5 If*/WN stars (e.g. Melnick 35) are very luminous, young 
stars with 150±30M G . However, some Of/WN stars are substantially lower in luminos- 
ity/mass (e.g. Sk-67° 22, Melnick 51), with correspondingly larger ages, even though 
these may be spectroscopically indistinguishable from other examples, as illustrated in 
Fig.EJ 

The incidence of 02-3.5 If*/WN stars in the LMC is significantly higher than 
in the Milky Way. Radiatively driven winds of Galactic stars would be expected to be 
modestly higher than LMC counterparts, so one would predict a slightly higher percent- 
age of 02-3.5 If*/WN stars in the LMC with respect to the Milky Way. In fact, 02- 
3.5 If*/WN stars comprise 7% of the 106 WN-type stars in the LMC dBrevsacher et all 
versu s only 2% of the highly incomplete 175 WN stars compiled by lvan der Huchtl 
lQl1 . l2006l) . Since transition spectral types arise preferentially in very massive stars, 
one would expect them predominantly in regions of the highest star formation. In- 
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Figure 3. Spectrograms of transition Of/WN stars dCrowther & Walbornll201 ll) 



deed, the 30 Dor regio n of the LMC dominates Of/WN statistics in the Local Group 
dCrowther & WalbornllloTTh . 



3. Death of Very Massive Stars 

A natural question relating to VMS is whether they would follow the usual path to core- 
collap se supernovae (CC SNe) or explode prematurely as pair-creation supernovae (PC- 
SNe)? iHeger et all (l2003h concluded that metal-free, single massive stars with 140-260 
Mq would explode as PCSNe. However, unbiased transient surveys have recently iden- 
tified exceptionally bright supernovae in the local universe, some of which have been 
attributed to PC SNe from initially ~ 200 M n stars with more modest metal-deficiencies 
(e.g. SN2007bi lGal-YametaDl2009h . 



VMS mode l s hav e been calculated throughout their post-main sequen ce evolution 
using l Vink et all ( 2001 ) mass-loss prescriptions for the main sequence and lNugis & Lamers 



42000D for the post -main sequence Wolf-Rayet phase, the results of which are pre- 
sented in Fig. 01 H-deficient CCSNe are predicted for 100-300 M Q stars at solar and 
LMC metallicities, whereas CO core masses of 60-130 M Q are obtained for rotating 
150-200 M Q stars at SM C metallicity. T herefore, VMS at low metallicity may indeed 
produce PCSNe. Indeed, Ouim by et al.1 d201 ll) have identified a class of lum inous H- 



deficient supernovae located in faint, metal-poor host galaxies. lOuimby et all attributed 
such bright SNe to a strong interaction between the CCSN of a very massive star and 
a H-free shell produced by violent pulsations, that was perhaps initiated by the pair 
instability. 
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Figure 4. Final CO mass versus initial mass for VMS at SMC, LMC and so- 
lar metallicities , with the domain of PCSNe shown in grey, b ased upon mass-loss 
prescriptions of IVink et al.l d2001l) and[Nugis & Lamersl (12000) for, respectively, the 
H-rich and H-poor phases (from Yusof et al. in prep.). 



However, these predictions are very sensitive to mass-loss prescriptions, especially 
for the post-main sequence phase. For the Wolf-Rayet phase separate expressions are 
adopted for WN and WC s tars, in which mass-loss rates are expressed in terms of lu- 
minosity and composition (|Nugis & Lamers 2000, eqns. 20-21). Such calibrations, 
based on results for Wolf-Rayet stars at solar compositio n, imply a factor of -30 in- 
crease in mass-loss rate from the H-rich dVink et al.ll200ll) to the H-dencient phase of a 
300 M Q star at SM C metallicity. To illus trate the sensitivity, let us alternatively adopt 
equation 25 from lNugis & Lamersl (120001) . albeit modified to allo w for the M oc Z" o 7 de- 
pendence of mas s-loss upon ambient ( Fe-peak) metallicity, Zp e (ICrowther et all 2002; 
Vink & de Koterll2005l : ICrowtheril2006l) . i.e. 



logMAMoyr" 1 ) = -5.7 + 0.88 log(M/M ) + 0.7 log(Z Fe /Z Fe , G ) 



This would exceed the I Vink et all d200ll) prediction by only a factor of 2 for the case of 



a SMC metallicity 300 M Q star, leading to significantly higher CO masses than those 
presented in Fig. 0] raising the possibility of PCSNe from VMS progenitors at higher 
metallicity. 
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Discussion 

Krysztof Stanek Is there any photometric variability information for your most mas- 
sive stars?: 

Paul Crowther Not that we know, although their location at the centres of very crowded 
clusters makes photometric studies challenging (R136al and R136a2 are separated by 
only 0.1 arcsec): 



